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Abstract 

The five different CP conserving amplitudes for the decays K — ► 3tt are calculated 
using Chiral Perturbation Theory. The calculation is made to next-to-leading order 
and includes full isospin breaking. The squared amplitudes are compared with the 
corresponding ones in the isospin limit to estimate the size of the isospin breaking 
effects. In this paper we add the radiative corrections to the earlier calculated m u — 
and local electromagnetic effects. We find corrections of order 5-10 percent. 



PACS numbers: 13.20.Eb; 12.39.Fe; 14.40.Aq; 11.30.Rd 



1 Introduction 



The non-perturbative nature of low-energy QCD calls for alternative methods of calculating 
processes including composite particles such as mesons and baryons. A method describing 
the interactions of the light pseudoscalar mesons (K, it, 77) is Chiral Perturbation Theory 
(ChPT). It was introduced by Weinberg, Gasser and Leutwyler [TJ Ej and it has been 
very successful. Pedagogical introductions to ChPT can be found in The theory was 
later extended to also cover the weak interactions of the pseudoscalars J5J, and the first 
calculation of a kaon decaying into pions (K — ► 2tt,3tt) appeared shortly thereafter [B]. 
Reviews of other applications of ChPT to nonleptonic weak interactions can be found in 



A recalculation in the isospin limit of K — > 2n to next-to-leading order was made in 
(El Ej and of K — > 3tt in j2l QDj • In |Hj also a full fit to all experimental data was made 
and it was found that the decay rates and linear slopes agreed well. However, a small 
discrepancy was found in the quadratic slopes and that is part of the motivation for this 
further investigation of the decay K — > 3tt in ChPT. 

The discrepancies found can have several different origins. It could be an experimental 
problem or it could have a theoretical origin. In the latter case the corrections to the 
amplitude calculated in [0] are threefold: strong isospin breaking, electromagnetic (EM) 
isospin breaking or higher order corrections. These effects have been studied in many 
papers for the K — > 2tc decays, references can be traced back from [TT]. For K —>■ 3n 
less work has been done. In the strong isospin and local electromagnetic corrections 
were investigated and it was found that the inclusion of those led to changes of a few 
percent in the amplitudes. The local electromagnetic part was also calculated in in 
full agreement with our result after corrections of some misprints in ^Hj • 

In this paper we add also the radiative corrections, i.e. the nonlocal electromagnetic 
isospin breaking. The full (first order) isospin breaking amplitude to next-to-leading chiral 
order is thus calculated, and we will try to estimate the effect of this in the amplitudes. 
A new full fit, including also new experimental data f3J E]> has to be done to answer 
the question whether isospin breaking removes the problem of fitting the quadratic slopes. 
This, together with a study of models for the higher order coefficients, we plan to do in an 
upcoming paper. 

Other recent results on K — > 3tt decays can be found in In [H] Nicola 

Cabibbo discusses the possibility of determining the ao — a% pion scattering length from 
the threshold effects of K + — ► 7r 7r 7r + . He gives an approximate theoretical result with 
very few unknown parameters. We have a possibly better theoretical description of these 
effects but it includes more unknown parameters. In [15] an attempt was made to calculate 
the virtual photon corrections to the K + — > 7r°7r°7r + decay. Our result disagrees with the 
result presented there. 

The outline of this paper is as follows. The next section describes isospin breaking in 
more detail. In section[H]the basis of ChPT, the Chiral Lagrangians, are discussed. Section 
0] specifies the decays and describes the relevant kinematics. The divergences appearing 
when including photons are discussed in section |3J In section |H1 the analytical results 
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are discussed, section [7| contains the numerical results and the last section contains the 
conclusions. 

2 Isospin Breaking 

Isospin symmetry is the symmetry under exchange of up- and down-quarks. Obviously 
this symmetry is only true in the approximation that m u = and electromagnetism is 
neglected, i.e. in the isospin limit. Calculations are often performed in the isospin limit 
since this is simpler and gives a good first estimate of the result. 

However, to get a precise result one has to include isospin breaking, i.e. the effects 
from m u ^ and electromagnetism. Effects coming from m u ^ we refer to as strong 
isospin breaking and include mixing between 7r° and rj. This mixing leads to changes in 
the formulas for both the physical masses of 7r° and rj as well as the amplitude for any 
process involving either of the two. For a detailed discussion see |17j . 

The other source is electromagnetic isospin breaking, coming from the fact that the up- 
and the down-quarks are charged, which implies different interactions with photons. This 
part can be further divided in local electromagnetic isospin breaking and explicit photon 
contributions (radiative corrections). The former are described by adding new Lagrangians 
at each order and the latter by introducing new diagrams including photons. 

Our first calculation of K —>■ 3tt || was done in the isospin limit. In the next paper, 
[T2] , we included strong and local isospin breaking (there collectively referred to as strong 
isospin breaking) and we now present the calculation including all isospin breaking effects. 

3 The ChPT Lagrangians 

The basis of our ChPT calculation is the various Chiral Lagrangians. They can be divided 
in different orders. The order parameters in the perturbation series are p and m, the 
momenta and mass of the pseudoscalars. Including isospin breaking also e, the electron 
charge, and the mass difference, m u — rrid, are used as order parameters. All of these 
are independent expansion parameters. We work to leading order in m u — and e 2 but 
next-to-leading order in p 2 and m? . For simplicity we call in the remainder terms of order 
p 2 , m 2 , e 2 and m u — leading order, and terms of order p 4 ,p 2 m 2 , m 4 , p 2 e 2 , m 2 e 2 , 
p 2 (m u — rrid) and m 2 (m u — ma) next-to- leading order. 

3.1 Leading Order 

The leading order Chiral Lagrangian is usually divided in three parts 

£-2 = C-S2 + C-W2 + C-E2, (1) 

where £32 refers to the strong AS = part, C W2 the weak AS = ±1 part, and Ce2 the 
strong-electromagnetic and weak-electromagnetic parts combined. For the strong part we 
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have 



C S 2 = ^(u^ + x + ) 



(2) 



Here (A) stands for the flavour trace of the matrix A, and Fq is the pion decay constant 
in the chiral limit. We define the matrices u^, u and x± as 



where the special unitary matrix U contains the Goldstone boson fields 



(3) 



U = exp i^y-M 



M 



( >S + > VT+ K + \ 



(4) 



The formalism we use is the external field method of [2], and to include photons we set 



X = 2£? 



V 



m„ 



\ 



and £> M 17 = <9 M £7 - ze QA^U - ie UQA 



(5) 



m s ) 

where A^ is the photon field and 




-1/3 J 

We diagonalize the quadratic terms in (j2J) by a rotation 

7T° = 7r 3 cos e + 7/s sin e 
77 = — 7T3 sin e + r/ 8 cos e . 

where the lowest order mixing angle e satisfies 

tan(2e) = V?>- 



(6) 



(7) 



m d - m u 



2 m s - m u - m d 
The weak part of the Lagrangian has the form ^H] 



C 



W2 



G 8 {A 32 u^) + G' 8 (A 32X+ ) + G 27 t ij ' kl (A ijU ,){A kl u») 



+ h.c. 



(9) 



The tensor t %h has as nonzero components 



21,13 _ ^13,21 



t 



+23,33 _ +33,23 



3 
1 



+22,23 _ ^23,22 



.23,11 _ +H,23 



1 

3 



(10) 



and the matrix A y - is defined as 

Ay = uXijU 1 , (Xij) ab = S ia 5 jb . (11) 
The coefficient C is defined such that in the chiral and large N c limits Gg = G27 = 1, 

C=-\ %V ud C = -1.06 • 1(T 6 GeV- 2 . (12) 

Finally, the remaining electromagnetic part, relevant for this calculation, looks like (see 

e.g. ng) 

C E2 = e 2 F*Z(Q L Q R ) + e 2 F^TQ R ) (13) 
where the weak-electromagnetic term is characterized by a constant Ge (g e wkGs in [T9]). 

T = G s F 2 A 32 + h.c. (14) 

and 

Ql = uQu\ Q r = u ] Qu. (15) 
3.2 Next-to-leading Order 

The fact that ChPT is a non-renormalizable theory means that new terms have to be added 
at each order to compensate for the loop-divergences. This means that the Lagrangians 
increase in size for every new order and the number of free parameters rises as well. At 
next-to-leading order the Lagrangian is split in four parts which, in obvious notation, are 

£4 = ^54 + C-Wi + C-S2E2 + C-W2E2{G%) . (16) 

Here the notation (G$) indicates that here only the dominant Gs-part is included in the 
Lagrangian and therefore in the calculation. 

These Lagrangians are quite large and we choose not to write them explicitly here since 
they can be found in many places E3 El EU 1221 DZ3 I2H] • For a list of all the pieces 
relevant for this specific calculation see ^2]- Note however that four terms producing 
photon interactions should be added to Cwa in [J2j. The two new terms in the octet part 
are 

n u i ( a 32 {/r, UflUv }) + n 15 i (a 32 vr«*> ( i7 ) 

and in the 27 part 

D 13 itv> kl (A ijUl t)(A kl K/f ]> +D 15 it i i> kl (A ij u l ,u„){A kl f!?), (18) 

where 

/£" = uF^v) + vlF^u, F» v = eQ{d»A v -d u A»). (19) 
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3.2.1 Ultraviolet Divergences 

The process K — > 3tt receives higher-order contributions from diagrams that contain loops. 
The study of these diagrams is complicated by the fact that they need to be defined pre- 
cisely. The loop-diagrams involve an integration over the undetermined loop-momentum 
q, and the integrals are divergent in the q — > oo ultraviolet region. These ultraviolet diver- 
gences are canceled by replacing the coefficients in the next-to-leading order Lagrangians 
by the renormalized coefficients and a subtraction part, see |HJE] and references therein. 
The divergences can be used as a check on the calculation and all our infinities (except the 
ones left since the G 2 r-pa,rt in £w2E2 is not known) cancel as they should. 

3.2.2 Loop Integrals 

The prescription we use for the loop integrals can be found in many places, e.g. [24J. The 
only one needed in addition to the ones given there is the one-loop three point function 

C(m 1 ,m 2 ,m 3 ,p 1 ,p 2 ,p 3 ) = - / , . , 7—5 27-77 ^ 2777 7> 2\ > (^) 

1 J (2ir) a (jr — mf) {{p — p\Y — mi) ({p — p?,y — m 3 ) 

where p 3 = pi+p 2 . For its numerical evaluation we use the program FF j2H]- This program 
also deals with possible infrared divergences consistently. 



4 Kinematics 

There are five different CP-conserving decays of the type K — >• 3n (K~ decays are not 
treated separately since they are counterparts to the K + decays): 



000] ) 



K L (k) - 7t ( Pl )n (p 2 )n (p 3 ), [A 

K L (k) - n + ( Pl )n-(p 2 )n°(p 3 ), [A L + _ ] , 

K s (k) - n+( Pl )n-(p 2 )n°(p 3 ), [A s + _ ] , 

K+(k) - 7r ( Pl )7r (p 2 )7r + (p3), [A)o+] , 

K+(k) - 7r + ( Pl )vr + (p 2 )7r-(p3), , (21) 

where we have indicated the four-momentum defined for each particle and the symbol used 
for the amplitude. 

The kinematics is treated using 

si — (k — pif , s 2 = (k-p 2 ) 2 , s 3 = (k-p 3 ) 2 . (22) 

The amplitudes are expanded in terms of the Dalitz plot variables x and y defined as 

S3 — So S 2 — S\ 1 . . , . 

y = — 2 — ' x = — 2 — ' s o = o ( Sl + 32 + S3 > ■ ^ 23 ) 

m% + m;+ 3 
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Figure 1: Bremsstrahlung, the emission of an extra final-state photon. 

The amplitude for Kl — » ttVtt is symmetric under the interchange of all three final state 
particles and the one for K$ — > 7r + 7r~7r° is antisymmetric under the interchange of 7r + and 
tt~ because of CP. The amplitudes for Kl — ► 7r + 7r~7r°, K + — > 7r + 7r + 7r~ and K + — > 7r 7r 7r + 
are symmetric under the interchange of the first two pions because of CP or Bose-symmetry. 



5 Infrared Divergences 

In addition to the ultraviolet divergences which are removed by renormalization, diagrams 
including photons in the loops contain infrared (IR) divergences. These infinities come from 
the q — > end of the loop-momentum integrals. They are canceled by including also the 
Bremsstrahlung diagram, where a real photon is radiated off one of the charged mesons, see 
Fig. It is only the sum of the virtual loop corrections and the real Bremsstrahlung which 
is physically significant and thus needs to be well defined. We regulate the IR divergence 
in both the virtual photon loops and the real emission with a photon mass m 7 and keep 
only the singular terms plus those that do not vanish in the limit m 7 — > 0. We include the 
real Bremsstrahlung for photon energies up to a cut-off u and treat it in the soft photon 
approximation. 

The exact form of the amplitude squared for the bremsstrahlung diagram depends on 
which specific amplitude that is being calculated. For K + {k) — ► 7r c '(pi)vr \p 2 )7i'^ '(P3) it can 
be written in the soft photon limit (see e.g. [2o] ) 



\A\ 2 
\ a \bs 



\ a \lo e 



d 3 q 1 



(2vr) 3 2q 



E 

A=0,1 



k ■ p 3 ■ 



q ■ k 



(24) 



where \A\lo is the lowest order isospin limit amplitude. The number of terms inside the 
parentheses is the number of charged particles in the process and the sign of those terms 
depends both on the charge of the radiating particle and on whether it is incoming or 



outgoing. Writing out the square and using 



A=0,1 



r(A) p (A) 



-9v», you get 



Ul 



BS 



\A\ 



LO 1 



d 3 q 1 



(2vr) 3 2q 



k' 2 



+ 



Pi 



2p 3 -k 



(q-k) 2 (q-p 3 ) 2 (q-k)(q-p 3 ) 



(25) 
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To solve the first integral term, place the vector k along the z-axis, i.e. 

k = (k°, 0, 0, k z ) and (k ■ q) 2 = {k°q° - k z q z ) 2 . 
Changing to polar coordinates that part of the integral now looks like 



- \A\ 



m 



K 



LO { 



8tt 2 



J dqd(cos9) 



(k°E 1 - k z q cos #) 2 ' 



(26) 



(27) 



where k 2 = m 2 K , q = E 1 and q z = qcos9 have been used. Solving the d(cos9) part is now 
straightforward and leads to 



- \A\ 



m 



LO< 



K 



8tt 2 



dq 



1 



k z \k°E 1 - k z q k Q E 1 + k z q J ' 



(28) 



Putting the two terms on a common denominator and changing variable to E~ leads to 



-1^1 



2 2 m K 

L ° 6 4tt 2 



dE. 



7 E2(P) 2 -(E2-m 2 )(^) 2 ' 



(29) 



where u is the photon energy above which the detector identifies it as a real external photon. 
We are only interested in the result in the limit m 7 — > 0, so it's enough to consider 



12 d E 



1 



\LO' 



4vr 2 



7 m 2 ,E 7 ' 



which gives the result 



\LO g7r2 



e 2 , uj 2 
log 



mi 



(30) 



(31) 



In a similar way one gets the result for the mixed term 



/ 



d 3 q 1 



(2tt) 3 2q 



P3 ■ k 



x s s 3 - m\ - ml 



UJ 



where 



(q-k)(q-p 3 )\ Ait 2 m K m^{\ - x 2 s ) 

- Am K m w / (s 3 - (m K - m n ) 2 ) - 1 

- 4m K m n / (s 3 - (m K - ^tt) 2 ) + 1 



log x s log — = Ii R (m K , m n , s 3 ) , 



(32) 
(33) 



In order to obtain the correct imaginary part we use the ie-prescription, which means 
S3 = s 3 + is. 

For the other amplitudes the calculations are similar and the resulting bremsstrahlung 
amplitudes are 

(34) 

,2 



I a l I 2 _ n 



\A\ 



-0\BS 



-0\LO 



4vr 2 



mi 



(35) 



«\/\/\# m #v/v/n# 



Figure 2: The tree level diagrams for K — > 3n. A filled square is a weak vertex, a filled 
circle a strong vertex, a straight line a pseudoscalar meson and a wiggly line a photon. 



\A 



S |2 
-0lB5 



I 4 I 2 



I /\S |2 J_ 



l^oo+lio^ 



log— T - 7/^(77^, 771^, S 3 ) 



mi 



CO 



log— 7 - I IR (m n} m K ,s 3 ) 



mi 



(36) 
(37) 



I A, 



|2 

IBS 



I A, 



lLO 



4vr 2 



2 log — - Im(m w , m K , s x ) - Ii R (m n , m K , s 2 ) 



mi 



+Im(ml, m 2 K , s 3 ) - Im(ml, m 2 n , si) - li R [m\, m 2 w , s 2 ) 
+/ /R (m^,m^,s 3 ) 



(3J 



When using the above, the divergences from the explicit photon loops cancel exactly. 

A similar problem shows up in the definition of the decay constants since we normalize 
the lowest order with F n + and F K +. Our prescription for the decay constants is described 
in App. |XJ 



6 Analytical Results 
6.1 Lowest order 

The four diagrams that could contribute to lowest order can be seen in Fig. El 

However, the two diagrams including photons turn out to give zero. This is obviously 
so for K + — > 7r + 7r°7r° and Kl — > 7r 7r 7r° since the 77r°7r° vertex vanishes as a consequence 
of charge conjugation. 

The reason why it vanishes for the other decays is somewhat more subtle and is the same 
as why the lowest order result for K -—>■ tt£ + £~ vanishes [2Zj- When doing a simultaneous 
diagonalization of the covariant kinetic and mass terms quadratic in the pseudoscalar fields, 
including those of the weak lagrangian Cy/2-, p 2 -terms of the form d^Kd^ir are absent 
and all weak vertices involve at least three pseudoscalar fields. This result should not 
change as compared to our calculation where the weak Lagrangian was not included in the 
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II — O ■ () 

— fX — fX 

Figure 3: Examples of diagrams of next-to- leading order with no photons. An open square 
is a vertex from Cy/i or Cw2E2, an open circle a vertex from C S4: or Cs2E2, a filled square a 
vertex from Cy/2 or Ce2 (AS 1 = 1) and a filled circle a vertex from C$2 or £#2 (AS* = 0). 

diagonalization. Thus in our case, the two diagrams on the right in Fig. El will together 
give zero contribution. 

This means that the lowest order result in the full isospin case is the same as when just 
including strong and local EM isospin breaking. This result we published before, the full 
expressions can be found in [12J. 



6.2 Next-to-leading order 

There are 51 additional diagrams contributing to next-to-leading order. They can be 
divided in three different classes and examples will be shown of each class. It should be 
noted that the argument in the previous subsection is not valid at this order. There now 
exist Kit^i vertices. The reason for this is that one can not diagonalize simultaneously all 
terms with two pseudoscalar fields when going to next-to-leading order. 

The first class of diagrams are the 13 which do not include explicit photons. They are 
the ones used in our earlier papers ^] and a complete list of them can be found there. 
Some examples are shown in Fig. El 

The second class of diagrams are the ones with a photon running in a loop. There are 
18 of these and some examples can be found in Fig. |3J Their evaluation is the main new 
result of this paper. They are also responsible for the infrared divergences discussed in 
Sect. The first diagram is an example where the photon is the only particle in the loop, 
a photon tadpole diagram. These vanish in dimensional regularization when only singular 
and nonzero terms in the limit m 7 — > are kept. 
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Figure 4: Examples of diagrams with photons in the loops. A filled square is a weak 
vertex, a filled circle a strong vertex, a straight line a pseudoscalar meson and a wiggly 
line a photon. 

The last class of diagrams is the ones with tree level photon propagators, 20 in total. 
They are photon reducible, i.e. if we cut the photon line the diagram falls apart. They 
are infrared finite and some examples can be seen in Fig. It turns out that for realistic 
values of the input parameters this class of diagrams give a negligible contribution to all 
K — > 3tt processes. 

We work to first order in isospin breaking, i.e. as soon as e 2 is present we set m u = m<f, 
m n + = m n o = m n and m K + = m K o = m K . Even then, the resulting amplitudes at next- 
to-leading order are rather long, and it does not seem very useful to present them here 
explicitly. However, the full expressions for the amplitudes are available on request from the 
authors or can be downloaded j2Hj- It should be noted that the amplitude for K L — > 7r°7r 7r° 
was given in |12j . It was stated that it was the full isospin breaking amplitude since no 
explicit photon diagrams can contribute to this process. This is not completely true, the 
amplitude will change indirectly through the definition of F n + and F K +, see App. 1X1 

The next-to-leading order amplitudes can in principle include all the low-energy coeffi- 
cients (LECs) from Cwa, £>S4, £-W2E2 and Csiei- The coefficients from the strong and elec- 
tromagnetic part of the Lagrangian are treated as input, which leaves G 8 ,G 2 7,N^, . . . , iVJ" 3 , 
D\,D T 2 ,D\, . . . , Dlj, D^q, . . . , D r zl and Z[, . . . , Z'[ 4 as undetermined. In total 41 unknown 
parameters. However, all of these do not appear independently, i.e. they multiply the same 
type of term, e.g. m\ or e^m^. It turns out, as discussed in [T2"] . that there are 30 inde- 
pendent combinations, denoted by K\ . . . K%q. For the 11 combinations already appearing 
in the isospin limit see and the 19 additional ones for the isospin breaking case can 
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XL SI 



Figure 5: Examples of diagrams with photon propagators. An open square is a vertex from 
Cwi or Cw2E2, an open circle a vertex from £54 or Cs2E2, a filled square a vertex from 
Cw2 or Ce2 (AS" = 1) and a filled circle a vertex from £g 2 or Ce2 (AS 1 = 0). A straight 
line is a pseudoscalar meson and a wiggly line a photon. 
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Table 1: The various input values used. 

be be found in [TS]. In addition to the 30 combinations found in [T2], four new coeffi- 
cients show up when including photons: N[ A) N"[ 5 , D[ 3 and D[ 5 . These all come from the 
third class of diagrams where a tree level photon is present. The four coupling constants 
show up in precisely the same combinations in K —>■ tt£ + £~ and can thus be determined 
experimentally in other decays. We therefore treat them as input. 

7 Numerical Results 

7.1 Experimental data and fit 

A full isospin limit fit was made in j2] taking into account all data published before May 
2002. One of the motivations for this continued investigation of isospin breaking effects is to 
see whether isospin violation can solve the discrepancies in the quadratic slope parameters 
found there. A new full fit will be done in an upcoming paper. The data from ISTRA+ 
[TB*] and KLOE |14j . which appeared after [5], will then also be taken into account. We do 
not present a new fit in this paper since estimates of the new combinations of constants 
should be done before attempting a full fit. 

7.2 Inputs 

The input values we use are presented in Table [TJ 
7.2.1 Strong and Electromagnetic Input 

There are different ways to treat the masses, especially in the isospin limit case. In jH] the 
masses used in the phase space were obtained from the physical masses occuring in the 
decays. However in the amplitudes the physical mass of the kaon involved in the process 
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was used and the pion mass was given by m 2 = § Hi=i 3 wi 2 ^ with % = 1,2,3 being the 
three pions participating in the reaction. This allowed for the correct kinematical relation 
S\ + $2 + S3 = m 2 K + 3m 2 to be satisfied while having the isospin limit in the amplitude but 
the physical masses in the phase space. The results in |Hj were obtained with the physical 
mass for the eta. Results with the Gell-Mann-Okubo (GMO) relation for the eta mass in 
the loops gave small changes within the general errors given in jU]. 

In the decays here, we work to first order in isospin breaking. We have rewritten explicit 
factors of m u — in terms of sin e according to 

2 

m u - m d = --^= (2m s - m u - m d ) sine. (39) 

In general we use the physical masses of pions and kaons in the loops but as soon as a 
factor of sine or e 2 is present we use a common kaon and a common pion mass. This 
simplifies the analytical formulas enormously. The kaon mass chosen is the mass from the 
kaon in the decay and the pion mass used is 3m 2 = J2i m ^ with i = 1,2,3 the three pions 
in the final state, i.e. the mass we used in the isospin limit case. For the eta mass we use 
in general the GMO mass in the loops but with isospin violation included, 

2 1 

m % = 3 ( m K+ +m 2 K0 -m 2 + ) + -mlo. (40) 

The possible lowest order contributions from the eta mass have been removed from the 
amplitudes using the corresponding next-to- leading order relation as described in |12j . 

The strong LECs L\ to L§ from Csa as well as sine come from the one-loop fit in [T7j . 

The constant Z from Ce2 we estimate via 

Z=^- 2 « + -<o), (41) 

which corresponds to the value in Tabled The higher order coefficients of Cea, K\ . . . Kn, 
are rather unknown. Some rough estimates exist but we put them to zero here. 

The IR divergences are cancelled by adding the soft-photon Bremsstrahlung. We have 
used a 10 MeV cut-off in energy for this and used the same cut-off in the definition of F v + 
and F K +. 

The subtraction scale fi is chosen to be 0.77 GeV. 
7.2.2 Weak Inputs 

The coefficients contributing in the isospin limit from Cw2 and Cwa are taken from the fit 
in jU]. The values of Gg, G27 in Tableware taken from 9J as well. We use as a reasonable 
estimate for Ge the value presented in |2"9j . 

No knowledge exists of the values of K u . . . iT 30 , so they are set equal to zero at \i = 
0.77 GeV. Tests were also made assigning order of magnitude estimates to them. This 
imparted changes in the isospin breaking corrections similar in size to those of the loop 
contributions. 
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7.2.3 Input relevant for the Photon Reducible Diagrams 

The two new constants -D13 and -D15 are set to zero since no knowledge exist of their values. 
Note that in order to contribute at all, they have to get values orders of magnitude larger 
than the expected size. 

The two other new constants, iV^ and -/V 15 , can be determined from K — > nl + l~ decays. 
For these decays, the branching ratios can be expressed as [2"T] 

BR(K + -> vr+e+e - ) = (3.15 - 21. lw+ + 36.1m; 2 ) • 1(T 8 \C G 8 /(9 ■ • GeV~ 2 ) | 2 
BR(K S -> 7r°e + e-) = (3.07 - 18. 7w s + 28Aw 2 s ) ■ 10" 10 |CG 8 /(9 • 10" 6 • GeV~ 2 )| 2 

(42) 

Using the measured central values |30] 

BR(K + -> vr+e+e-) = 2.88 ■ 10~ 7 , BR(iT 5 -> TrVe") = 3.0 ■ 1(T 9 (43) 

one gets the results 

w ( + ] = 1.69, = -1.10 and = 1.93, w^ 2) = -1.28. (44) 

These constants, w+ and ws, can then be written in terms of both strong and weak 
low-energy constants [2*7j . 

w + = |(4vr) 2 [iV[ 4 ( /U )-iV[ 5 ( /U ) + 3^( /U )]-iln^^ 
3 3 n A 

ws = ^(4vr) 2 [2iV[ 4 ( /U ) + iV[ 5 ( /U )]-iln^f, (45) 

and from the knowledge of £g(/x) |31j . one gets the values listed in Tabled for one choice 
of signs in eq. (jUj). 

Note that the influence of this class of diagrams is not visible in the figures nor in the 
results shown in the tables for any of the possible signs chosen in eq. (JUJ). The magnitude 
of these constants needs to be increased significantly in order to have a visible impact on 
our numerical results. 

7.3 Results with and without isospin breaking 

The results we will present here is a comparison between the squared amplitudes and the 
decay rates in the isospin limit and including first order isospin breaking. The full squared 
amplitudes over the decay region is a 3-D plot with the two different cases plotted over 
phasespace. This is however very difficult to read, and instead we will present comparisons 
along three slices of these 3-D plots as explained below. We also present the corrections 
for the Dalitz plot parameters. 

In Table 121 we present the values of the amplitudes squared in the center of the Dalitz 
plot, i.e. for x = y = 0. We show the results in the isospin limit from [Hj, with the strong 



14 





Centralvalue 


Iso [9 


Strong (T21 


Full 


K L -> ttVtt 
Ki — > 7r + 7r~7r° 
— > 7r + 7r~7r° 

-»■ 7T 7r 7T + 
— > 7T + 7r + 7T- 


6.74 • 10" 12 
7.46 • 10~ 13 


0.93- 10~ 13 
3.72 • 10~ 12 


6.97- 10" 12 
7.66 • 1(T 13 


1.01 ■ 10~ 12 
4.00 ■ 10~ 12 


7.04- 10~ 12 
7.88 • 10- 13 


1.03 ■ 10~ 12 
4.14- 10~ 12 



Table 2: Comparison of the central values of the amplitudes squared in the isospin conserv- 
ing case (Iso), including strong and local electromagnetic (Strong) and full (Full) isospin 
breaking. 

and local electromagnetic isospin breaking included [T5] and with full isospin breaking 
included. 

Similarly, in Table El we present the integrated decay rates in the isospin conserving 
case jHj, the one with strong and local electromagnetic isospin breaking included ^2] and 
with all isospin breaking effects included. There are here in principle problems with an 
infinite correction when a charged two pion system is at rest. The effects of the electro- 
magnetic interaction can then become very large from terms containing logarithms of the 
pion velocity. This is where the Coulomb interaction dominates and it should then really 
be resummed to all orders. In order to avoid this problem we have introduced the cut-off 
Eq. It means that we only integrate the phase space over the part where 

y/Ti > 2m, + E c . (46) 

Due to the way we have chosen the pion mass, the Coulomb problem only shows up for 
the decay K + — > 7r + 7r + 7r~ where the systems of two charged pions can be at rest or at 
very low relative velocity at the edges of phase space. The places where this happens are 
indicated by the large dots in Fig. El The choice of the pion masses in K^ s ~^ 7r + 7r~7r is 
such that the Coulomb threshold is slightly outside the physical phasespace. It turns out 
that in this case the part of the correction that includes the Coulomb singularity is rather 
small. We have therefore not included any corrections for it in the results presented. 

In Fig. El we also show the phase space boundaries for the five different decays and the 
three curves along which we will show results for the squared amplitudes with and without 
isospin breaking. The three curves are x — 0, y — and x = V3y. In Fig. [7] to Fig. [TT] we 
then plot the five different squared amplitudes along these curves as a function of r, where 
r = ±^Jy 2 + y an d the sign is chosen according to 

!y, x = 

x/V3, y = (47) 
yy/2, x = yV3. 

Note that for all but A+_ the squared amplitudes are normalized to their value at the 
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Decay Rate 








Iso 




r\t TTlTl (X 
kJ IjI Wilt! 


Full 


hiQ [Mevj 


Kt - 


~ /I /I /I 


2.65-10 


-18 


2.74 • 10" 18 


2.77-10' 


-18 




K L - 




1.63-10- 


-18 


1.67- 10- 18 


1.72 • 10- 


-18 




K s - 




3.1 • 10- 


21 


3.2 • 10- 21 


3.3 • 10" 


21 




K+ - 


-> 7T 7r 7T + 


9.11 ■ 10" 


-19 


9.84- 10" 19 


i.oo-io- 


-18 




K+ - 




2.97-10" 


-18 


3.19 ■ 10- 18 













2.95-10" 


-18 


3.17- 10- 18 


3.28-10- 


-18 


1 






2.91 • 10" 


-18 


3.13 • 10- 18 


3.24- 10- 


-18 


2 






2.72 • 10' 


-18 


2.93 • 1Q- 18 


3.03-10 


-18 


5 



Table 3: Comparison of the decay rates in the isospin conserving case (Iso), including 
strong and local electromagnetic (Strong) and full (Full) isospin breaking. The Coulomb 
cut-off used, Ec, is explained in the text. 

center of the Dalitz plot. A comparison of the central values themselves is shown in Table 121 

We also calculate the changes in the Dalitz plot distribution parameters. These are 
defined by 

2 

= 1 + gy + hy 2 + kx 2 . (48) 

The isospin breaking corrections to these parameters are given in Table |3J The amplitude 
for K$ — > 7T + 7T-7r is parametrized via 

^+-o = isx - £ s xy ■ (49) 

We now discuss the results in somewhat more detail. In general the results are of a size 
as can be expected from this type of isospin breaking. They are of order a few, up to 11% 
in the amplitudes squared outside the Coulomb region. The isospin breaking corrections 
tend to increase all decay rates somewhat and this will in a fit be compensated by small 
changes in the values of the Ki compared to the results of [9 . The number of significant 
digits quoted in Table [21 is higher than the expected precision of our results, but the trend 
and the general size of the change compared to the isospin conserving results are stable 
with respect to variations in dealing with the eta mass (physical or GMO). 

For Kl — > 7r 7r°7r° the central value of the amplitude squared increases by about 4.5%. 
In this case we have because of the symmetry of the final state that g = and k = h/3. 
The quadratic slope decreases by about 5% but the total variation over the Dalitz plot is 
small so the total decay rate increases by about 4.5% as well. This decay is the one which 
has most variation in the amplitude when changing how one deals with the eta mass. The 
extreme case we have found was that this effect completely cancelled the change from 
isospin violation, but the relative change due to isospin breaking remained similar. Note 



A(s 1 ,s 2 ,s 3 ) 
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Figure 6: The phase space boundaries for the five different decays and the three curves 
along which we will compare the squared amplitudes. The points where two charged pions 
have low relative velocity for K + — > 7t + 7t + tt~ are indicated by the large dots. 



Decay 


Quantity 


Iso [9J 


Strong 


Full 




-> 7T U 7T U 7r U 


h 


-0.0072 


-0.0068 


-0.0068 


K L 


— > 7T + 7T~7T° 


9 


0.673 


0.683 


0.677 






h 


0.085 


0.089 


0.088 






k 


0.0055 


0.0057 


0.0057 


K s 


— > 7T + 7r _ 7T 


Is 


3.4 ■ 10~ 8 


3.4 • 10~ 8 


3.5- 10~ 8 






& 


-0.2 ■ 10" 8 


-0.2 ■ 10~ 8 


-0.2 • 10~ 8 


K + 


-> 7T 7r 7T + 


9 


0.635 


0.619 


0.619 






h 


0.074 


0.071 


0.071 


K + 


— > 7r + 7r + 7r _ 


9 


-0.215 


-0.211 


-0.201 






h 


0.012 


0.012 


0.008 






k 


-0.0052 


-0.0050 


-0.0037 



Table 4: Comparison of the Dalitz plot distribution parameters in the isospin conserv- 
ing case (Iso), including strong and local electromagnetic (Strong) and full (Full) isospin 
breaking. 
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Figure 7: Comparison of Ki — > 7r°7r 7r° with and without isospin breaking. 

the scale in Fig. when viewing the result. The changes compared to ^2] are entirely due 
to the photon loop corrections to F n + and F K +. 

In Fig. [7| one can also clearly see the thresholds induced by the difference between 
and m^ introduced when isospin invariance is broken. These thresholds correspond 
to a new process being allowed where two of the neutral pions are produced through an 
intermediate on shell state with one positive and one negative pion. 

The squared amplitude K L — > 7r + 7r~7r° increases by about 5.5% with very little variation 
with the eta mass treatment. The decay rate increases by the same amount. The changes 
in the Dalitz plot slopes are rather small as can be judged from Fig. |H| The marginal 
differences compared to g quoted in [5j are due to a slightly different fitting procedure to 
the amplitudes squared. 

For the decay Kg — > 7r + 7r~7r° the amplitude in the center of the Dalitz plot vanishes 
because of the symmetries. The amplitude and the slopes increase by about 3% as can be 
seen in Fig. |U] and the tables. 

The decay K + — > 7r°7r°7r + has a large increase. The squared amplitude in the center 
changes by about 11%. The linear slopes decrease somewhat leading to an increase of 
about 10% to the total decay rate when compared with the isospin conserved case. This 
is shown in Fig. El 

The amplitude for K + — > 7r 7r°7r + is also calculated in |IBJ. Our results don't agree 
with the numerics presented there. We find an increase in the amplitude while there a 
decrease is found. In [TH] a different choice of lowest order was made than here and in 
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Figure 8: Comparison of Kl — > 7r + 7r tt with and without isospin breaking. 
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Figure 9: Comparison of K s — > n + n tt° with and without isospin breaking. 
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Figure 10: Comparison of K + — > 7r°7r°7r + with and without isospin breaking. 

After taking that difference into account, we still disagree significantly with the 
numerical results of [To] . 

The decay K + — > 7r + 7r + 7r~ has a change of about 11% upwards in the center of the 
Dalitz plot. The slopes decrease somewhat. The decay rate can only be compared when 
the Coulomb region is excluded from the comparison but the total change is also about 
11%. 

The conclusions above do not change qualitatively when we give the K\ a value of about 
0.001 and the new isospin breaking Ki a value relative to G§ and G27 of 0.01. However the 
changes induced by these values are numerically significant. They can be of the order of 
10%, largest for K L -> 7r°7r 7r°. 

It should be noted that the mentioned changes are with the values of Ki determined 
from the isospin conserving fit in [0J. A new determination including isospin breaking 
effects is planned in an upcoming paper. 

8 Conclusions 

We have calculated the K — > 3ir amplitudes to next-to-leading order (p 4 ,p 2 m 2 , m 4 ,p 2 e 2 ) in 
Chiral Perturbation Theory. A similar calculation was done in [§| in the isospin limit, and 
in ^2] including strong isospin breaking, but we have now included full isospin breaking. 
The motivations for this are both because it is interesting in general to see the importance 
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Figure 11: Comparison of K + — ► 7r + 7r + 7r with and without isospin breaking. 



of isospin breaking in this process, but also to investigate whether isospin violation will 
improve the fit to experimental data. Discrepancies between data and the quadratic slopes 
from ChPT were found in |9] , and isospin breaking may be the cause of this. 

We have estimated the effects of the isospin breaking by comparing the squared am- 
plitudes with and without isospin violation. The effect seems to be at 5-10% percent level 
in the amplitudes squared. To investigate if this removes the discrepancies found in [U] a 
new full fit has to be done, also including the new data [131 E] published after 0. This 
is work in progress and will be presented in the future paper Isospin Breaking in K — > 3tt 
Decays III. 
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A The Decay constants F n + and F K +. 

We have chosen to normalize our lowest order contribution with Fq /(F^+F k +). F t + and 
F K + are the pion and kaon decay constants respectively. Including isospin breaking they 
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Figure 12: Diagrams for the decay constants. An open square is a vertex from Cwa or 
£-W2E2, a filled square a vertex from Cw2 or Ce2 (AS 1 = 1) and a filled circle a vertex from 
C$2 or Ce2 (AS = 0). A straight line is a pseudoscalar meson and a wiggly line a photon. 



are determined from the diagrams in Fig. ED and the resulting expressions are 

1 



F 1 + 



F 2 



F 2 e 2 (4/3 AT[ + 4/3 fT 2 r + 10/9 K\ + 10/9 ^ + 2 K[ 2 



<9 <9 

1 o "757 2 2 2\ 2 15 1 2 2 2\ n 15 ( '2- 2 2\ 2 



" 7 > 
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+ 1/2 A(rr4 + ) + 1/4 A(m 2 K0 ) + A(m 2 ) (3/8 - 3/4 -j= 
These formulas agree with known results jSl EH] • 



(A.2) 
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The above formulas are infrared divergent when m 7 — > 0, and the standard way to deal 
with this is the same as in the amplitudes, i.e. adding a bremsstrahlung diagram. We have 
chosen to just add a term including a cut-off scale for the bremsstrahlung photon, 

-|^Uo^log4> ( A -3) 

which cancels the dependence on the photon mass and therefore removes the divergence. 
The scale uif is set to 10 MeV. 
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